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As developing national and international space ventures , which seek to employ 
NASA's Deep Space Network (DSN) for tracking and data acquisition , evolve , it 
is essentia/ for navigation and tracking system analysts to evaluate the operational 
capability of Deep Space Station antennas. To commission the DSN for use in 
tracking a highly eccentric Earth orbiter could quite possibly yield the greatest 
challenges in terms of slewing capability; certainly more so than with a c/eep-space 
probe. This article focuses on the determination of the maximum slew rates needed 
to track a specific high Earth orbiter , namely, the Japanese MUSES-B spacecraft 
of the Very Long Baseline Interferometry Space Observatory Program. The results 
suggest that DSN 34-m antennas are capable of meeting the slew rate requirements 
for the nominal MUSES-B orbital geometries currently being considered. 


I. Introduction 

The Institute for Space and Astronautical Science 
(ISAS) has proposed the implementation of an orbiting ra- 
dio telescope observatory mission known as the Very Long 
Baseline Interferometry (VLBI) Space Observatory Pro- 
gram (VSOP) for launch in the mid-1990s timeframe. The 
program is structured to be an international collaboration 
of the VLBI radio astronomy communities of the United 
States, Europe, and Australia, as well as other interested 
nations, with the intention of providing a worldwide ob- 
servatory network. 

ISAS and NASA will undertake a cooperative effort 
for precise tracking of VSOP’s MUSES-B spacecraft by 
commissioning NASA’s Deep Space Network (DSN) to- 
gether with various Japanese ground stations. 


Due to the highly eccentric orbit proposed for 
MUSES-B, concern arose as to whether or not DSN 
34-m antennas would be capable of slewing fast enough 
to follow the spacecraft. An effort w r as undertaken to for- 
mulate an analytic model for the determination of antenna 
slew angles (i.e., azimuth and elevation angles), as well as 
antenna slew rates, and produce a computer program for 
generating their corresponding time histories. Analysis of 
different scenarios together with knowledge of the opera- 
tional slewing capability of the Deep Space Station (DSS) 
antennas would thereby yield tracking feasibility results. 

This article details the modeling and analysis pro- 
cess of this effort as well as results and implications of a 
few selected mission scenarios. Results indicate that DSN 
34-m antennas are indeed capable of meeting slewing 


13 


T 


requirements for MUSES-B’s nominal mission orbit and, 
more importantly, further preliminary studies can now be 
readily performed for future planned missions, 

II. Modeling Assumptions and Analysis 

In order to arrive at an analytical formulation of the 
problem, the spacecraft orbit was modeled with two-body 
dynamics. The Earth’s precession of equinoxes was not 
modeled. Classical orbital elements were assumed time- 
invariant as were the cylindrical coordinates of the station 
position. 1 The derivation of the azimuth and elevation 
angles in a (north, east, down) topocentric-horizon refer- 
ence frame yielded inverse trigonometric functions of the 
station-to-spacecraft slant range vector. The actual equa- 
tions are omitted here for brevity. Once the analytical 
expressions for the slant range vector were derived, ana- 
lytical expressions for the time-rates-of-change of azimuth 
and elevation could be determined. It should be noted 
that expressions for the slant range vector components 
were derived as functions of the time-varying parameters 
of station right- ascension, 9(t) y and spacecraft eccentric 
anomaly, E(t ); together with the time-invariant parame- 
ters of the classical elements, station position components, 
and associated direction-cosine terms. 

Despite the fact that a simple model was used for this 
study, the analytical expressions for azimuth and elevation 
angle determination and their corresponding time-rates-of- 
change were too lengthy to obtain a time history of results 
by hand. The expressions were mechanized in a computer 
program to facilitate rapid evaluation of different cases. 
Input/Output specifications are listed below. 

A. Inputs 

The inputs are r,, X 3ta , z 0 , 9 g0 , w®, a, e, f), i, ui, p®, 
E 0 , t 0 , dt, tsrop where 

r s — station distance from Earth’s spin axis (km) 
A 5ia = station longitude (east of prime meridian) 
(deg) 

Zq == s_t at ion distance above equatorial plane (krn) 
0 gQ = Greenwich mean sidereal time (sec) 

= Earth’s rotation rate (rad/sec) 
a — semi-major axis (km) 
e = eccentricity 

Q = longitude of the ascending node (deg) 

1 T. Moyer, “Station Location Sets Referred to the Radio Frame,” 

JPL IOM 314.5-1334 (internal document), Jet Propulsion Labo- 

ratory, Pasadena, California, February 24, 1989. 


i — inclination (deg) 
w = argument of perigee (deg) 

/i® = Earth’s gravitational parameter (km 3 /sec 2 ) 
E 0 = initial spacecraft eccentric anomaly (sec) 
to — specified epoch time (sec) 
dt — time increments (sec) 
tsTOP = program stop time (sec) 

B. Outputs 

The outputs are t y 6 — 9(t), E = E(t) } Az = 
Az(0, Ej;p), El = El(9, E y t;p), t, 0 = 6{t ), E = E(t), 

Az = Az(9,E,t;p), El = El(6,E,t,p) ( Az ) max , (£/) max 
where 

t — time (hrs) ^ - - 

9 ~ station right-ascension (deg) 

E = spacecraft eccentric anomaly (deg) 

Az = azimuth angle (deg) 

El = elevation angle (deg) 

Az = azimuth rate (deg/sec) 

El = elevation rate (deg/sec) 
p — input parameter set 

III. Results and Implications 

A case-by-case synopsis is provided for three relevant 
mission scenarios. Tables of input parameters are included 
as are plotted slew angle time histories. The referenced 
epoch was taken to be January 1, 1995 at o hr 0 mm 0 sec and 
duration of each run was seventy- two hours. The orbit 
period in all cases was approximately 6.06 hours. 

A. Case 1: (MUSES-B Nominal Orbit; DSS 10) 2 

The Deep Space Station represented in this case is 
located at Goldstone, California. Input parameters are 
listed in Table 1 while plotted results of slew angles are 
illustrated in Figs. 1 and 2. 

Maximum slew rates for this case were found to be 
approximately 

0.14 deg/s in azimuth 
0.10 deg/s in elevation 

B. Case 2: (MUSES-B Nominal Orbit; DSS 40) 

This scenario was identical to Case 1 with the excep- 
tion that the Deep Space Station represented is located at 

2 A. Konopliv, “Preliminary Orbit Determination Analysis for the 
VSOP Mission,” JPL IOM 314.5-648 (internal document), Jet 
Propulsion Laboratory, Pasadena, California, February 9, 1989. 
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Canberra, Australia. The input parameters are listed in 
Table 2. Plotted results of slew angles are illustrated in 
Figs. 3 and 4. 

Maximum slew rates for this case were found to be 
approximately 

0.14 deg/s in azimuth 
0.09 deg/s in elevation 

C. Case 3: (MUSES-B Nominal Orbit; DSS 60) 

This scenario was also identical to Case 1 with the 
exception that the Deep Space Station represented is lo- 
cated at Madrid, Spain. The input parameters are listed 
in Table 3. Plotted results of slew angles are illustrated in 
Figs. 5 and 6. 

Maximum slew rates for this case were found to be 
approximately 


0.14 deg/s in azimuth 
0.08 deg/s in elevation 

IV. Conclusions 

Maximum slew rates for a DSN 34- m antenna are 
roughly 0.40 deg/s in both azimuth and elevation and , 
hence, the 34-m antenna would certainly be capable of 
meeting slewing requirements for any of the case scenarios 
described above. 

Station view periods can also be determined by exam- 
ining elevation angle time histories that are generated by 
the computer program. Fig. 7 provides an example of this 
capability, which represents Cases 1 through 3 employing 
an elevation cutoff angle of ten degrees above the local 
horizon. 
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Table 1. Input parameter set for Case 1 


7*3 (km) 

^ sta (deg) 

-0 (km) 

5204.00 

243.1105 

3677.05 


^so(deg) 

(rad/sec) 

- 


- • • 


100.175 

7.292116 X 10- 5 





a (km) 

e 

Q(deg) 




16878.00 

0.5629 

0.00 




'(deg) 

w(deg) 

" - • r - ' 

iLl jr.-'- .11 



46.00 

0.00 


Table 3. Input parameter set for Case 3 

(km 3 /sec 2 ) 

Eo (deg) 





3.9860045 X 10* 

0.00 


r 3 (km) 

Aita(deg) 

^o(km) 




4862.45 

355.752 

4115.11 

<o(sec) 

dt(sec) 

tsrop{ sec) 




0.00 

1200.00 

2.592 x 10 5 

5 so(deg) 

(rad/sec) 






100.175 

7.292116 X 10~ 5 





a(km) 

e 

fi(deg) 




16878.00 

0.5629 

0.00 




'(deg) 

u/(deg) 





46.00 

0.00 


Table 2. Input parameter set for Case 2 

M®(km 3 /sec 2 ) 

£J 0 (deg) 





3.9860045 X 10 5 

0.00 


r*3 (km) 

Aata(deg) 

zq (km) 




5205.25 

148.9813 

-3674.75 

fo(sec) 

dt( sec) 

*STOp(sec) 




0.00 

1200.00 

2.592 x 10 5 

^jo(deg) 

u/0(rad/sec) 






100.175 

7.292116 x 10- 5 


a(km) 

e 

Q(deg) 

16878.00 

0.5629 

0.00 

2 (deg) 

w(deg) 


46.00 

0.00 


^9 (km 3 /sec 2 ) 

^o(deg) 


3.9860045 x 10 s 

0.00 


io(sec) 

dt( sec) 

f STOP (sec) 

0 . 0 Q 

1200.00 

2.592 x 10 5 



ELEVATION, deg AZIMUTH, deg 


400 



Fig. 1. Azimuth versus time (Case 1). 



Fig. 2. Elevation versus time (Case 1). 
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ELEVATION, deg AZIMUTH, deg 



Fig. 5. Azimuth versus time (Case 3). 



Fig. 6. Elevation versus time (Case 3). 
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Fig. 7. DSS view periods with 10-deg lower elevation cutoff, 
epoch 1995-01-01 00:00:00.0000. 
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